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ABSTRACT 

Segue 1 is the current best candidate for a “first galaxy”, a system which experienced only a 
single short burst of star formation and has since remained unchanged. Here we present possible star 
formation scenarios which can explain its unique metallicity distribution. While the majority of stars 
in all other ultra-faint dwarfs (UFDs) are within 0.5 dex of the mean [Fe/H] for the galaxy, 5 of the 7 
stars in Segue I have a spread of A[Fe/H] > 0.8 dex. We show that this distribution of metallicities 
canot be explained by a gradual build-up of stars, but instead requires clustered star formation. 
Chemical tagging allows the separate unresolved delta functions in abundance space to be associated 
with discrete events in space and time. This provides an opportunity to put the enrichment events 
into a time sequence and unravel the history of the system. We investigate two possible scenarios for 
the star formation history of Segue 1 using Fyris Alpha simulations of gas in a 10^ Mq dark matter 
halo. The lack of stars with intermediate metallicities —3 < [Fe/H] < —2 can be explained either by a 
pause in star formation caused by supernova feedback, or by the spread of metallicities resulting from 
one or two supernovae in a low-mass dark matter halo. Either possibility can reproduce the metallicity 
distribution function (MDF), as well as the other observed elemental abundances. The unusual MDF 
and the low luminosity of Segue 1 can be explained by it being a first galaxy that originated with 
Mvir 10^ Mq at z ~ 10. 

Subject headings: early universe; galaxies: dwarf; galaxies: abundances; galaxies: formation; stars: 
abundances; dark ages, reionization, first stars 


1. INTRODUCTION 

The formation and evolution of the first stars and 
galaxies can be studied through observations of nearby 
dwarf galaxies. Stars preserve the chemical signatures 
at the time of their formation, such that galaxies with 
an early truncation of their star formation history can 
be used to investigate the conditions of the high-redshift 
universe. These systems are called “first galaxies” or 


“fossil galaxies” (Ricotti & Gnedin 

2005 

Bovill (fe Ri- 

cotti 

120091 IBromm & Yoshida||2011 

F're 

3el & Bromm 

2012 

Frebel et al. 2014 Bland-Haw 

thorn et al. 2015), 


frozen in time for > 10 Gyr. The best candidates for such 
systems are the ultra-faint dwarfs (UFDs, L < 10® Lq) 
discovered in the past decade. Most UFDs contain only 
old, metal-poor s tars and appear to have s imple star for¬ 
mation histories (Brown et al.|20l'2 20141, making them 


excellent candidates tor probing the chemical signatures 
of the first generations of stars. 

There is not yet a consensus on the definition of a “first 
galaxy”, but most authors exclude the ~ 10® Mq mini¬ 
halos in which the first Population HI stars formed and 
require a first galaxy to produce and retain a long-lived 


stellar system (Bromm & Yoshida 2011 Bland-Hawthorn 


I d.webster@physics.usyd.edu.au 


et al. 20151. Frebel & Bromm (2012) defined a first 
galaxy as one which had only a single short burst of star 
formation. Such a system would be enriched only by Pop¬ 
ulation HI stars and would therefore lack the signatures 
of Type la and AGB star enrichment events, resulting 
in enhanced alpha element abundances and low or zero 
neutron-capture element abundances. 

The best candid ate to date for a first galaxy is Segue 1 
dFrebel et al . 120141) . Segue 1 was discovered b y|Belokurov| 
etaC [20071 ! and confirmed as a galaxy by jCeha et al. 


(120091), w ho showed that it wa s dark matter-dominated 

who showed that it had 
The mean metallicity of 


as well as Norris et al. (2010b) 
a large spread in meta 
Segue 1, [Fe/H] 


Jlicity. 

—2.7 made it the most metal-poor 
known galaxy, although the recently di scovered Retic¬ 
ulum 2 has a similar me an metallicity (Koposov et al. 
2015b| Shimo n et ar||2015 ). High resolution spectroscopy 


(Frebel et al.|2014|) showed that the 7 brightest red giants 
in Segue 1, with mietallicities ranging from [Fe/H] = —3.8 


] 


to —1.4, showed [a/Fe] « 0.5 and suppressed [Sr/H] and 
[Ba/H], suggesting a maximum of one r-process or weak 
s-process event. This implies that no stars formed from 
gas enriched by AGB stars, suggesting that star forma¬ 
tion in Segue 1 lasted less than the lifetime of 7-8 Mq 
stars, ~ 30 — 50 Myr. 

The metallicity distribution function (MDF) of Segue 1 
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is also unusual. While all other known UFDs have metal- 
licity distributions that peak near the mean metallicity 
and have only a tail of metal-poor or metal-rich stars, 
Segue 1 shows three stars with [Fe/H] « —3.6, two with 
[Fe/H] « —2.4 and two with [Fe/H] « —1.5. Five of 
the seven stars are more metal-poor or metal-ri ch than 


90% of obse rved stars in the 6 UFDs observed by Brown 
et al.| (2014). In this work, we seek to explain the unusual 


MUF of Segue 1 by exploring possible scenarios for its 
star formation history. This can give us an insight into 
the conditions of star formation in the early universe, 
as well as suggesting new avenues for determining which 
systems are likely to be Hrst galaxies. 

In the lo cal universe, star fo rmation is observed only 
in clusters ( Lada fc Lada||20d3 |. There is some evidence 
to suggest that t his should also be true for the high- 


redshift universe (Clark et al. 

T 


_ ,,2008 Larsen et al. 2012, 

Karlsson et al.||20l^ 2013). Clusters have homogeneous 


abundances as long as they are not self-enriching, mean¬ 
ing that the timescale for its star formation must be less 
than the time of the first supernova. This is suppor ted 


Silva et al.[[2006 

2007 

Feng & Krumholz 

2014). 

'i'he chemicai homogeneity ot clusters 

has led to the 


the search for stars with similar abundances in a range of 
elements, suggesting that they originated from the same 
cluster. Separate unresolved delta functions in chemical 
abundance space are associated with events which are 
discrete in space and time. This provides an opportunity 
to put these events into a time sequence and determine 
t he detailed enrichment h i story o f a system. 


Bland-Hawthorn et al. (2010a I tailored che mical tag¬ 
ging to dwa rf galaxies, and it was applied by |Karlssoii| 


et al. (2012) to Sextans, the best candidate for a cluster 


signature m a dwarf galaxy. They found that three stars 
with [Fe/H] = —2.7 contained remarkably similar Mg, 
Ti, Cr and Ba abundances. If correct, this group of stars 
is more metal-poor than any other known star cluster in 
the Milky Way. 

Recently, 17 ultra-faint satellites have been discov¬ 


ered in the Southern sky (Koposov et al. 2015a Drlica- 


Wagner et al. 20151. It is not yet conhrmed whether they 


are all Uh'Us, as it is possible that some are globular clus¬ 
ters. The only systems to date with [Fe/H] data from 
spectroscopy are Reticulum 2 and Horologium I, which 
are confirmed as UFDs from measurement s of the veloc¬ 
ity dispersion and the spread in [Fe/H] ([Walker et al. 


2015 Simon et al. 2015 Koposov et al. 2bl5bp . Retic¬ 


ulum 2 is approximately tied with Segue 1 as the most 
metal-poor galaxy, with [Fe/H] = —2.7. Reticulum 2 is 
one of the 6 objects described by [Koposov et al. (2015a I 
as having sizes and luminosities similar to Segue 1. How- 
ever, as we will see in Section]^ the MDF of Reticulum 2 
is more similar to the more luminous UFDs than to that 
of Segue 1. 

In Section [^ we discuss the MDFs produced under 
various assumptions for the star formation history of a 
UFD. We compare these to the MDFs of Segue I, Retic¬ 
ulum 2, and the entire UFD population. In Section [^ 
we present possible star formation histories for Segue 1, 
guided by Fyris alpha hydrodynamical simulations of 
small (Mvir = 10^ M©) dark matter halos. We discuss 


our results in Section |4| and summarise in Section [5] 

2 . THEORETICAL MDF 
2.1. Method 

In this section, we show MDFs resulting from various 
star formation histories and compare them to the MDF of 
Segue 1. Given that Segue 1 is believed to have formed 
stars for less than 50 Myr, it should have a relatively 
uncomplicated star formation history. 

The first two panels of Figure [^ correspond to models 
of continuous enrichment. In the hrst panel (labelled 
“Type H”) an equal amount of star formation occurs 
between each of 14 Type 11 supernovae. Between time t 
and t + At, 100 stars form with [Fe/H] for a given star s 
given by: 


[Fe/H](s) = log(nFe(s)/nH) - log(nFe, 0 /nH. 0 ) 


= log(- 


».«(.)) _,, ^ ( 1 ) 


HFeMii 


where Mh = 10 ® Mq and Mp ^., the initial mass of iron 
in the gas, is set such that [Fe/H] = —3.5 at t = 0, ^Fe 
is the atomic mass of iron, N is the number of supernova 
explosions per unit time, ype = 0.08 Mq is the iron yield 
of each (Type H) sup ernova, log(nFe, 0 /fiH, 0 ) = —4.55 
(Asplund et al. 2005) and X ~ A/'(0,0.3^) is a ran- 
dom number selected from a Gaussian distribution with 
y, = 0, a = 0.3, representing the approximate di spersion 
in metallicities i n the interstellar medium (ISM) (Feng & 
Krumholz 2014). NAt is set to equal 1, such that there is 


one supernova explosion between t and t + At, the period 
during which 100 stars form. The “Type H-|-Ia” model 
in the second panel is the same, except that for a single 
value of t, taken to be half the total time, ype = 0.8 Mq, 
corresponding to the yield of a Type la sup e rnova. The 
supernova yields are from Iwamoto et al. (1999) and 
Nomoto et al. (2006). The model was run 100 times 
and the average taken. 

The neutral hydrogen mass Mh = 10® Mq assumed is 
of the order of the amount of gas in the inner 100 pc of an 
Mvir = 10® ®“^ Mq dark matter halo with a baryon frac - 
tion of 10%. As shown in Bland-Hawthorn et al. (2015|), 
such a system can retain the majority of its gas and met- 
als in the face of a single supernova explosion. The gas 
initial constant metallicity [Fe/H] = —3.5, corresponds 
to the expe cted metallicity for such a system enriche d by 
a first star (Greif et al. |2007l [Frebel & Bromm|2'M^. At 
each timestep the h'e yield of a single Type II supernova 
is added to the total Fe of the gas, with the neutral hy¬ 
drogen content assumed to remain constant, such that 
the metallicity of the gas gradually increases. 

The second class of scenarios, shown in Panels 3-6 
(“One Burst”, “Two Bursts”, “Two Bursts 75/25” and 
“Three Bursts”) of Figure are simply sums of Gaus¬ 
sian distributions. In all cases, the scatter of metallici¬ 
ties within a burst is given by the standard deviation of 
the Gaussians, set to be 0.3 dex as for the continuous 
models. 

This gives a total of 6 scenarios: 


1. “Type II”: Continuous enrichment from Type H 
supernovae as in Equation I. 
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Fig. 1.— Theoretical normalised MDFs (filled histograms) resulting from a range of possible simple star formation histories, along with 
the normalised MDFs observed in UFDs (unfilled histograms). In the top six panels, the mean [Fe/H] ^ —2.7, as in Segu e 1 and Reticulum 
2. The lower three panels show the observed MDFs of in dividual UFDs, as well as the entire oliserved UFD population (|Kirby et al.|2013| 
[Brown et al.|2014| [Frebel et al.|20T4l [Simon et al.|2015| |. 


2. “Type II + Type la”: Continuous enrichment from 
Type II and Type la supernovae as in Equation I. 

3. “One Burst”: A single burst of stars with mean 
[Fe/H] = —2.67, represented by a single Gaussian. 

4. “Two Bursts”: The sum of two equally-weighted 
Gaussian distributions with means [Fe/H] = —3.17 
and —2.17, representing the stars being evenly di¬ 
vided between the bursts. 

5. “Two Bursts 75/25”: The sum of two Gaussians, 
one with mean [Fe/H] = —2.87 and one with —2.17. 
The first Gaussian is given a relative weight of 3 
times the second, such that this is a model in which 
75% of the stars form in the lower metallicity burst. 

6. “Three Bursts”: Three evenly-weighted Gaussians 
with mean [Fe/H] = —3.7, —2.5 and —1.6, repre¬ 
sentative of three bursts with the stars evenly di¬ 
vided between the bursts. 

The first two scenarios correspond to closed box mod¬ 
els, with no inflows or outflows, where the iron is mixed 
instantaneously with the neutral hydrogen gas. Star for¬ 
mation is terminated when [Fe/H] in the stars formed is 
approximately —2.7, which is the mean [Fe/H] observed 
in Segue 1 and Reticulum 2. This number of supernovae 


is consistent with the 150 0 Mf^ of star formati on im¬ 
plied for S egue 1 under a Kroupa (20011 IMF ( Frebel] 


et al. 2014). Assuming a specihc st^ formation rate of 

^ 2008 | 

ill take 


10 ^M[/yr ^kpc ^ (Bigiel et al. 


20141 as in the Garina dwarf, this wil 


Webster et al. 
100-200 Myr. 

The number of stars in a given [Fe/H] bin increases ex¬ 
ponentially with increasing [Fe/H] until the mean metal¬ 
licity at which star formation is switched off, because 
[Fe/H] is a logarithmic scale, and the Fe mass is assumed 
to increase linearly, with the hydrogen mass remaining 
constant. There are a few stars with higher metallicities, 
which result from the assumed dispersion of [Fe/H]. 

While Segue 1 shows no evidence of Type la enrich¬ 
ment, we include the second panel (“Type H -b Type 
la”) to show how Type la supernovae alter an MDF. 
The build-up of metallicities remains exponential for the 
low-metallicity stars which form first. However, a Type 
la supernova injects a large amount of Fe, resulting in a 
rapid increase in Fe which flattens the peak of the [Fe/H] 
distribution. 

Scenarios 3-6 represent distinct bursts of star forma¬ 
tion, with no self-enrichment within each burst and the 
gas considered to be well-mixed between each burst. 
They give an indication of the variety of qualitatively 
different MDFs which can be constructed using s imple 
bursts. This is motivated by Brown et al. (2014), who 
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used isochrone fitting to determine the star formation 
history of 6 UFDs, finding that they were best fit by two 
bursts of star formation, although in 3 of the galaxies, 
> 95% of the stars were in one of the bursts. Giving the 
bursts a non-zero duration did not improve the fit. Sce- 


nario 3 (“One Burst”) has been proposed as a scenario for 

a first galaxy ( 

Frebel & Bromm 

2012 

), while the fourth 

scenario ( Two Bursts”) is simi 

Lar to the 

Brown et al. 


Burits 75/25”) shows two bursts w ith the first contain- 


ing 75% of the stars, similar to the Brown et al. (20141 
star formation histories for Hercules and Leo IV. 

The lower three panels of Figure [l] show data from ob- 


se rvations of Segue 1 ( [Frebel et al.p014[ ) and Reticulum 
2 (Simon et al.|2015), as well as the combined data from 
the six U FUs observed by Brown et al. (2014) along with 
Segue 2 ([Kirby et al.||20l3|). Segue 1 and Reticulum 2. 
While Reticulum 2 could be explained by a number of 
possible scenarios, Segue 1 is highly unusual. The clos¬ 
est match would be a three-burst scenario with smaller 
dispersion than shown in the “Three Bursts” panel. The 
overall UFD population shares qualitative features with 
the continuous Type II -I- la model, but this could also 
be caused by it being a composite of a large number of 
bursts across many systems. 

2.2. Clustering in Segue 1 

Figuref^shows how unusual Segue 1 is compared to the 
overall IJFD population. For all UFDs except Segue 1, 
the majority of stars have metallicities within 0.5 dex of 
the mean. In Segue 1, 5 of the 7 stars have metallici¬ 
ties outside the range of > 90% of the stars in UFDs. 
As well as the extreme metallicities, any scenario to ex¬ 
plain the star formation history of Segue I needs to ex¬ 
plain the gaps in metallicity between the three stars with 
[Fe/H] « —3.6, the two with [Fe/H] « —2.4 and the 
two with [Fe/H] « —1.5. No such gap is seen anywhere 
else in MDFs of UFDs, with the exception of Canes Ve- 
natici H, which contains a single high-metallicity star. 
The fa voured star formati on history for Canes Venatici 
H from Brown et al. (2014) is two bursts, with the second 
burst separated troiii the first by 3.2 Gyr containing only 
5% of the stars, so this higher metallicity star is likely 
part of a small second burst. 

As a statistical test to compare a three cluste r model to 
an unclustered model, we use the Beale (1969) pseudo-F 
statistic: 


j.. J? - Jj (N - .2) 

H (AT-1) - (A^-C2)c”^/^ 

where J/ = - [lY and - (J-kY 

correspond to the squared errors of the unclustered and 
clustered models, are the cluster centroids, N is the 
number of data points, C2 is the number of clusters and 
p is the dimensionality of the data. The first term is 
a comparison of how well the data matches the models, 
while the second term includes only the number of clus¬ 
ters, number of observations and the dimensionality of 
the data, such that models with more parameters are 
penalised. 

For a three-cluster model with centroids at 
(—3.65,-2.36,-1.55) compared to a model with¬ 



[Fe/H] 


Fig. 2.— The normalised MDF of Segue 1 l|Frebel et al.||2014|l 
(blue) and the combined normalised MDF ot Bootes 1, Uoma 
B erenices. Canes Ve natici II, Herc ules. Leo IV and U rsa Major 

1 llBrown et al.|2014 ||. Reticulum 2 (|Simon et al.|2015|| and Segue 

2 jKirby et aL|2Ui^ (green). 


out clusters, F* = 274 « F(2,4)p=o.oooij showing that 
the three-cluster model is a better fit to the data. 
The three-cluster model also outperforms the best 
two-cluster model, with F* = 7 « F(l, 4)p=o.05- 

Gaps in the MDF can be produced by star formation in 
discrete bursts. Single-age bursts should produce groups 
of stars that form at about the same time and therefore 
have similar metallicities. However, even three bursts 
centered on the metallicity of the three groups in Segue 1 
will not necessarily produce an MDF like in Segue 1. 
While there may be gaps in the MDF, typically one or 
more of the groups will contain either no stars or only one 
star, which does not conclusively show that the groups 
exist. Put differently, it is possible that there are systems 
with a similar star formation history to Segue 1 where we 
will not observe one or more of the groups and the system 
will therefore not be recognised as having grouped star 
formation. 

There are good physical and observational reasons to 
expect that star formation should be clustered in time 
and/or place in low mass systems. In the l ocal uni¬ 


verse, stars are observed to form in clusters (Lada & 


Lada 2003). If this is also the case at the very low gas 


masses an d star forma t ion ra tes of the UFDs, as sug¬ 
gested by Clark et al. (2008), we would expect to ob¬ 
serve groups ot stars with the same chemical composi¬ 
tion. In larger systems, this effect is washed out because 
eventually stars exist at all metallicities. However, in a 
“first galaxy” such as Segue 1, there are no later gener¬ 
ations, meaning that the signature of the first c l usters 
are observed. As shown in Bland-Hawthorn et al. (2015) 
and Webster et al. (2014), in a small enough dark mat- 
ter halo, a single supernova can affect the entire system, 
with star formation only possible during sufficiently long 
breaks between supernovae. If there are a large num¬ 
ber of supernovae during a short period of time, the gas 
may be blown out completely, permanently ending star 
formation in the system. 

The unusual MDF of Segue 1 is likely connected to its 
status as a first galaxy, with both features caused by its 
inability to retain gas for more than ^ 50 Myr. This is 
plausible if Segue 1 had a mass Mvir ~ 10^ Mq when it 
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formed stars at z > 1 0, which would imply a current mass 
of Mvir = 10®“® Mq ( Webster et al.|[20T5 ). For Segue 1, 
the observed half-light mass Mhaif = 2 zb 1 x 10® Mf^ 
and half-light radius rhaif ~ 30 pc (Collins et al. 20141 
give a best fit mass assuming an Mh'W prohle iWvir ~ 
3 X 10® Mq, but with an uncertainty of more than an 
order of magnitude. 

Given that continuous enrichment is unlikely, Segue 1 
must have formed its stars in separate bursts. However, 
these bursts can not have been similar to those deter¬ 


mined by Brown et al. 

a range of | Fe/H| and Ict/Fe] an d likely had some self¬ 
enrichment (Webster et al.|20l5|. For Segue 1, the three 


(2014), which produced stars with 
^e] 


bursts each contain red giants which are consistent within 
the uncertainties with the stars within each cluster hav¬ 
ing the same [Fe/H]. 

2.3. Carbon 

In the Milky Way halo and in UFDs, a large propor¬ 
tion of very metal-poor stars are enhanced in carbon rel¬ 
ative to iron. This is seen in Segue 1, where 4 of the 
7 stars with measured abundances are carbon-enhanced 
metal-poor (CEMP) stars, with [C/Fe] > 0.7, although 
one of the four ap pears to have receiv ed carbon from a 
binary companion (Frebel et al. 2014). This star is ex¬ 
cluded from the discussion below. I'he scatter A [C/Fe] is 
> 2 dex, significantly larger than for heavier elements in 
Segue 1, suggesting that carbon enrichment was decou- 
pled from the enrichment of heavier elements (see also 


Frebel et al.||2014| |Ritter et al.||2015| [Sluder et al.||2015|) . 

I'igure f^shows the relationship between [C/HJ and 
[Fe/H], along with lines representing constant [C/Fe]. 
There is no clear trend of [C/H] with [Fe/Hj. This is sug¬ 
gestive that the stars with higher [Fe/H] in Segue 1 need 
not be a consequence of a large number of enrichment 
events, but are instead the result of inhomogeneous mix¬ 
ing. With few events, the inherent scatter that results 
from a single supernova becomes more important, as we 
will see in Section |3l 

Figure also provides some additional evidence for 
clustered star formation in Segue 1 in that despite the 
high overall scatter in carbon abundances, the two stars 
with [Fe/H] « —2.4 have similar [C/H], as do 2 of the 
3 stars with [Fe/H] « —3.5. The only outlier as classi¬ 
fied by the Fe-clusters is the star with [C/H] = —1.25 
at [Fe/H] = —3.57. This star is from a different obser¬ 
vational sample to the other s and is also an outlier in 
[Mg/Fe], [A l/Fe] and [Ni/Fe] ( Norris et al.||2010b{ Frebel 
et al.||2014|). 


Hitter et al. (2015) explain stars with enhanced [C/Fe] 
as resulting from differences in the entropy of the inner 
and outer regions of the supernova ejecta. In their sim¬ 
ulations, the gas which first returns to the center is defi¬ 
cient in the innermost 10% of the ejecta, which contains 
iron peak and a elements. The outer ejecta cools first 
and does not rise to the same radii as the inner ejecta. 
The outer ejecta collapses to the center of the system, 
resulting in carbon-enhanced, iron-poor star formation. 
In Segue 1, there is low variance in [X/Fe] abundances of 
stars for a and Fe-peak elements, but [C/Fe] does vary 
sig nihcantly between c lusters. This is consistent with 
the Ritter et al. (2015) single and clustered supernovae 
simulations. 



Fig. 3. — [C/H] vs [Fe/H] for Segue 1, with dashed lines denoting 
lines of constant [C/Fe], The empty circle indicates the star which 
has been carbon-enhanced by its binary companion, while the filled 
in circles are stars which likely retai n the abundance o f the gas in 
which they fo rmed. The data is fromlFrebel et al.|(|2014|l andlNorrisI 

[7ni] ( i2oiob| i. '-' '—' 


3. SCENARIOS FOR THE STAR FORMATION 
HISTORY OF SEGUE 1 

In the discussion above, we considered only the gaps 
in the MDF and were not concerned with the precise 
metallicities. We have not yet provided a solution as to 
how stars with metallicities as high as [Fe/H] ^ —1.5 
can form. No other UFD contains more than one ob¬ 
served star with [Fe/H] more than 1 dex higher than its 
mean [Fe/H]. This is despite the fact that unlike Segue 1, 
most of these systems likely had a longer star formation 
history, as they show evidence of enrichment from AGB 
stars and Type la supernovae. 

Reaching an average [Fe/H] « —1.5 with only Type H 
supernovae from the [Fe/H] « —3.6 in the lowest metal- 
licity three stars would require enrichment by « 70 events 
10® MfTi o f gas. The stellar mass of 1500 M^ in 


per 

Segue 1 given a Kroupa (|2001|) IMF therefor e implies too 
few e nrichment events. A top-heavy IMF (Geha et al. 
2013) could cause sufficient enrichment events, but it 
is likely that 70 supernovae per 10® Mq of gas within 
50 Myr would blow out most of the metals and much o f 
the gas ( Mac Low fc Ferrai^|1999 Webster et al.||20i4 ). 

Inhomogeneous mixing must therefore be at least part 
of the explanation for the highest metallicity stars. How¬ 
ever, gas with metallicity [Fe/H] = —1.5 rarely shows 
densities sufficient for star formation in g as accreted onto 
a dark matter halo dCreif et al.||2010t, or after a sin¬ 
gle supernova event ( Bland-Hawthorn et al.||2015 ). Any 
proposed star formation history must not only explain 
why there are no stars with metallicities between [Fe/H] 
= —1.7 and —2.3, or between [Fe/H] = —2.5 and —3.4, 
but also how stars with such a high metallicity can form 
despite the short timescale of star formation in the sys¬ 
tem. 

We consider three possibilities to explain the MDF. In 
the first, the clusters of higher metallicity stars form at 
the interface of two colliding supernova remnants, where 
dense gas is swept up and partially mixed with the ejecta. 
The second scenario has the entire system enriched to 
[Fe/H] « —2.5 by ~ 10 SNe, with no star formation 
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Fig. 4. — Two possible star formation histories for Segue 1. Both models start with a cluster or clusters of stars forming at [Fe/H] 5^ —3.5 
(green stars) and finish with stars with the metallicities observed in Segue 1 today. The upper track has a cluster of stars forming at 
[Fe/H] Pi —1.5 (red stars) at the collision of two supernova remnants. Subsequent supernovae enrich the gas to [Fe/H] Pi —2.5, at which 
time another cluster of stars forms (yellow stars). Star formation is then terminated within 50 Myr of the first cluster forming. In the lower 
track, the gas is first enriched to [Fe/H] —2.5 by ~ 10 supernovae, with little star formation between them. There is then a gap between 
supernovae in which gas returns to the center. The first gas to return is more metal-rich and forms a cluster of stars with [Fe/H] Pi —1.5 (red 
stars), with the cluster at [Fe/H] Pi —2.5 (yellow stars) forming later when the gas is more well-mixed. Star formation is then terminated 
either by supernovae or a global event such as reionization. 


possible in between, because each supernova temporarily 
pushes out the dense gas. The [Fe/H] « —1.5 stars form 
from the first dense gas that returns to the center in a 
pause between supernovae, with the [Fe/H] « —2.5 stars 
forming shortly after. Finally, we consider the possibil¬ 
ity of enrichment by exotic types of SNe such as pair- 
instability supernovae. 

3.1. Simulations 


Our hydrodynamical simulations are described in de- 


tail in 

Bland-Hawthorn et al. 

20151). We used the Fvris 

Alpha hydrodynamical code 

Sutherland 

2010), which 


solves ideal Euler hydrodynamical hows of gas with cool¬ 
ing and an equation of state appropriate to astrophysical 
gases in a fixed gravitational potential. A nested three- 
level grid was set up, with each level having 216^ cells 
and being a factor of 3 smaller than its outer containing 
level. The largest level covers the halo out to beyond 
the virial radius of 630 pc, while the level containing 
the inner region has a resolution of 1.4 pc per cell. At 
t = 0, a 25 Mq star is assumed to form in a low mass 
(Mvir = 10® — 10^ M©) dark matter halo. 

The effect on the gas, which was initially set up as a 


fractal interstellar medium in approximate equilibrium, 
was traced for 25 Myr with a resolution of 0.5 Myr. The 
models include the energy output of a 25 Mq star both 
during the main sequence lifetime and from the super¬ 
nova after 6 Myr. The ionisation phase from the precur¬ 
sor star was modelled using the MAPPINGS IV ionisa¬ 


tion code (Allen et al. 


tion structures, the A 


2008) for the thermal and ionisa- 


LAS9 atmospheric ctM ( Castelli 


& Kurucz 20041, and the Meynet & Maeder ] 2002p evolu¬ 


tionary tracks. T he full details are in Section 2 of[Web-| 


ster et al.| (20141. The supernova was modelled by in- 
serting a bubble of hot gas with an equivalent energy 
of 10®i erg. As the high-pressure bubble expands, its 
thermal energy is converted to kinetic energy and radia¬ 
tive losses. We assume the star explodes as a typical 


core-collapse supernova. Faint supernovae (e.g. Cooke & 


Madau 2014) and other exotic supernovae will be dis¬ 
cussed in Section 3.4. 

The ionisation and supernova from a 25 Mq star per¬ 
manently switched off star formation for systems with 
M < 10® ® Mq, but above this limit, dark matter ha¬ 
los retain sufficient baryons to continue forming stars. 
A star positioned away from the center couples less ef- 
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ficiently with the dense gas, with much of the energy 
escaping in the directions away from the center. This 
results in more gas retention and less enrichment. In a 
10^ Mq dark matter halo with a baryon fraction of 10%, 
star formation is switched off for ^ 10 Myr after the ex¬ 
plosion of a high mass star near the center, but continues 
almost unaffecte d if the progenitor st ar is > 80 pc away 
from the center ( Webster et al.|[20l4 ). 

In this section we consider two variants of the above 
model as approximations to two possible star formation 
histories of Segue 1. The first involves the collision of 
a centered and off-centered explosion, while the second 
models a single supernova after the gas has been enriched 
to [Fe/H] = -2.5. 


3.2. Merging supernova remnants 

One possible scenario for the star formation history 
of Segue 1 involves the highest metallicity stars form¬ 
ing at the site where two supernova remnants collide. A 
single supernova exp anding in an HII region can also in- 


singie supernova exp anding m an Mil region can also in¬ 
duce star formation (Whalen et al.|2008 |Nagakura et al. 
20091, but here we investigate the colliding case, which 


20091, but here we investigate the colliding 
allows for higher metallicities to be reached. This possi- 
bility was b r iefly c onsidered in relation to first galaxies by 


Greif et al. (20071, who noted that if a collision between 
SJN remnants occurred early in their evolution in an over- 
dense region, the density could become high enough for 
gravitational frag mentation and therefore st ar formation. 


As predicted by Ostriker & McKee (1988), a supernova 
remnant develops a dense shelf of gas when its age is ap¬ 
proximately equal to the cooling time of the shock-heated 


1 Rosen & Bregman 1995 

Korpi et al. 1999 

de Avillez 

2000 Joung & Mac Low 

200t)|) have found that when 


these shells collide, dense and cold filamentary clouds are 
formed due to thermal instabilities and supersonic turbu¬ 
lence, an effect first no ted for the case of two int eracting 
supernova remnants by Bodenheimer et al. (1984). In our 
simulations, given gas densities p > 1 cm“'^, which are 
observed in the swept-up gas, the cooling time is shorter 
than our 500 kyr time resolution and star formation is 
therefore possible. If sufficient metal mixing occurs in 
the swept-up gas, the stars will form at a high metallic¬ 
ity. To investigate this scenario, we use d the Fyris Alpha 
si mulations discussed briefly in Section |3.1| and in detail 
in|Bland-Hawthorn et al. (2015). 

'I'he scenario is illustrated in the upper path of Fig¬ 
ure 1^ At t = 0, there is a dark matter halo containing 
gas enriched to [Fe/H] = —4. This starting metallicity 
is uncertain and could be as high as [Fe/H] = —3.5. In 
this gas, a cluster of stars forms. This cluster includes 
stars with masses M « 0.75 Mq observed as red giants 
today, denoted by the green stars, along with at least 
one star with M > 8 Mq (black stars in Figure]^. At 
a similar time, a cluster forms in lower density gas well 
away from the center with at least one massive star. The 
two clusters contain the 3 stars with metallicities of ap¬ 
proximately [Fe/H] = —3.5 observed as red giants today. 

One supernova then occurs in each of the clusters de¬ 
scribed above. The inner supernova ejects 0.08 Mq of 
Fe into a region with 2 x 10^ Mq of gas, enriching it 
to [Fe/H] = —2.5. At the same time, the supernova 
further out, which occurs in less dense gas, has enriched 
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Fig. 5.— A central slice of the metallicity and density of the gas 
just before the collision of supernova remnants from central and 
off-centered explosions. The two remnants have swept up dense 
neutral gas. 


a region with 2 x 10^ Mq of gas to [Fe/H] = —1.5. 
These supernova remnants expand, sweep up gas and 
collide, inducing star formation in the swept-up gas. A 
cluster of stars forms with [Fe/H] = —1.5 (red stars in 
Figure]^. The density and metallicity in the simulations 
at the time at which the two supernovae remnants col¬ 
lide is shown in Figure [^ Note that these simulations 
are unable to describe the complex physics of colliding 
supernova remnants, so here we are interested only in 
establishing whether the scenario is plausible. 

Subsequent supernovae then enrich the gas to [Fe/H] = 
—2.5 and a cluster forms at [Fe/H] « —2.5 containing two 
low-mass stars (yellow in Figure® that will be red giants 
today. These stars must form less than 50 Myr after the 
first cluster of stars, because stars in Segue 1 do not show 
evidence of AGB or Type la enrichment. Finally, gas is 
blown out, terminating star formation. The gas can be 
blown out either by multiple closely-spaced supernovae, 
or as a result of the epoch o f reion ization. As noted by 


Nichols & Bland-Hawthorn (2011), heated gas can be 
removed through tidal stripping by the host galaxy, or 
ram pressure stripping by a tenuous a hot halo. 

In the merging supernova remnants scenario, the rea¬ 
son that stars form at [Fe/H] « —1.5 is that this is the 
metallicity of dense gas swept up and mixed with the 
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supernova ejecta. In the future we will seek to per¬ 
form hydrodynamical simulations of the effects of col¬ 
liding supernova remnants, including whether the con¬ 
ditions for star formation exist, but this is beyond the 
scope of the present work. Here we simply note that at 
the boundary of the remnants, gas with [Fe/H] similar 
to the highest metallicity stars can build up, such that 
given the right conditions, stars will form. The stars at 
[Fe/H] « —2.4 may be formed slightly later at another 
site along the boundary between the two remnants, or 
by the first dense gas that returns to the center of the 
system after ~ 10 Myr, which is more metal-rich than 
the mean metallicity, as we will see in Section [3.3| 

We therefore suggest the following explanation for the 
MDF in Segue 1. The first stars formed with [Fe/H] « 
—3.6, near the mean metallicity of the system when it 
formed. The first two SNe occurred at approximately 
the same time, separated by at least ^ 75 pc. Dense gas 
was swept up and partially mixed with the supernova 
ejecta. A shock formed and the gas cooled, resulting in a 
cluster of stars forming with [Fe/H] much higher than the 
mean metallicity of the system. All the star formation in 
the system could be complete within 10-15 Myr, before 
supernova feedback blew all the gas out. A top - heavy 
initial mass function as suggested by Geha et al. (2013|) 
is expected to pro duce 250-400 supernovae in Segue 1 
( Frebel et al.||2014 ), which could easily blow the gas out 
of an Mvir ~ 10 ^ Mq dark matter halo. 

This scenario is particularly interesting in light of the 


Ritter et al. (2015) simulations discussed in Section 2.4. 
They found enhanced carbon relative to iron for gas 
which falls into the center. The iron-rich ejecta has more 
entropy and reaches higher radii than the carbon-rich 
ejecta, resulting in the stars which form at the center be¬ 
ing carbon enhanced. Here we instead have stars forming 
at higher radii where two supernova remnants merge. If 
the iron has overtaken the carbon by the time of the colli- 

Fe/H] and low 
Fe/H] = -1.5 


sion, the stars which form will have high 
[C/Fe]. This is observed in the star at 
with C/Fe] ~ 0. It should be noted that [[Ritter et al. 
(2015) model lower-mass (10® M©) systems over longer 


time periods, but it is suggestive of a possibility for in¬ 
homogeneous enrichment. 


3.3. Two bursts 

Another possible star formation history for Segue 1 is 
shown in the lower track of Figure and can be sum¬ 
marised as follows: 

1. A cluster of second-generation stars formed at t = 0 
with [Fe/H] « —3.5 (shown in green in Figure]^. 

2. Within the 30-50 Myr allowed for star formation, « 
10 SNe resulted from the star formation of this first 
cluster, enriching the medium to [Fe/H] « —2.5. 
This supernova rate is too high for the gas to re¬ 
cover between events, so little or no star formation 
occurs during this time. 

3. Due to the stochasticity of supernova events at low 
star formation rates, there was a sufficiently long 
gap r ^ 10 Myr between supernovae, allowing the 
gas to return to the center. The first gas to reach 
the center can have metallicities as high as [Fe/H] « 


— 1.5 and reach sufficient densities that stars may 
form (red stars in Figure® Gas in different parts of 
the system remained at [Fe/H] « — 2.5 and formed 
another cluster (yellow stars in Figure]^. 

4. As a result of the star formation from these clus¬ 
ters, or due to the epoch of reionization or other en¬ 
vironmental effects, the star forming gas was evap¬ 
orated or blown out of the system, resulting in the 
termination of star formation. 


We investigated this possibility using the F yris Al¬ 
pha simulations describ ed in Section 3.1 and in jBlanJ] 
Hawthorn et al. ( 2015| ). However, instead of a uniform 


starting metallicity [h'e/H] = —4, we assumed a starting 
metallicity of [Fe/H] = —2.5, such that the simulation 
models only the final supernova in step 2 of the scenario 
discussed above. In the initial state of the simulation, 
we assume that all the gas has this starting metallicity. 
In reality, it is likely that the gas will not be this well- 
mixed, and that regions of higher metallicity will exist. 
These pre-existing inhomogeneities may assist in reach¬ 
ing metallicities higher than we find. 

We assume that a supernova occurs in a medium al¬ 
ready somewhat disrupted by t he pr evious supernovae. 
While Bland-Hawthorn et al. (2015) models a 25 Mq 
star, here the star will be lower mass, because ^ 10 su¬ 
pernovae from the cluster have already occurred, sug¬ 
gesting a time scale of 20 — 30 Myr. This implies a mass 
of < 15 M 0 . However, the effect of the previous super¬ 
novae will be such that the surrounding medium is likely 
to be largely evacuated. As an approximation, the super¬ 
nova is therefore assumed to occur in the same medium 
as that for a 25 Mq star. 

The effect of the supernova is shown in Figure]^ The 
gas that returns to the center after ~ 10 Myr is en¬ 
riched to higher [Fe/H] then the overall average. While 
the metallicity of the gas in the central 20 pc has mean 
[Fe/H] « —2.0, there are pockets of dense gas with metal¬ 
licities [Fe/H] > —1.6. For example, 11 Myr after the su¬ 
pernova, there is a clump of gas with a mass of 200 Mq 
and [Fe/H] = —1.6 within an approximately spherical 
region with a radius of 7 pc. 

The third cluster observed today may form at a sim¬ 
ilar time in one of the lower metallicity clumps. The 
combined effect of the supernovae from these two newly- 
formed clusters, or the onset of the epoch of reionisation, 
may then eject all the gas in the system such that the 
three clusters described above contain all the stars that 
ever formed in Segue 1. 

This scenario provides a natural explanation of the ob¬ 
served MDF. The first cluster at [Fe/H] « —3.5 is re¬ 
sponsible for the three lowest metallicity stars observed 
as red giants today. The supernovae resulting from this 
cluster enrich the ISM to [Fe/H] « —2.5, while blowing 
the gas out of the center such that stars with metallic¬ 
ities in the range —3.5 < [Fe/H] < —2.5 are prevented 
from forming. The other two groups of stars form dur¬ 
ing a gap between supernovae, before star formation is 
permanently switched off. 

The abundances of elements relative to iron should 
then be close to the average Type H supernova yields, 
although there may be a bias towards the higher mass 
stars, which explode first. The first 10 supernovae may 


































9 



-200 -150 -100 -50 0 50 100 150 200 

x(pc) 



-200 -150 -100 -50 0 50 100 150 200 


x(pc) 


Fig. 6.— A central slice of the metallicity and density of the gas 
11 Myr after a supernova, the time at which the first dense gas is 
returning to the center. 

not be the entire supernova output of the first cluster, in 
which case the lower mass supernova progenitors would 
not yet have exploded by the time the higher metallic¬ 
ity clusters in Segue 1 form. This is indeed seen, with 
[a/Fe] - 0.5. 

Under this scenario, part of the explanation for Segue 1 
is that it has just the right mass and supernova rate 
that star formation can be temporarily switched off by 
supernova feed back without the gas being blown away 
completely. Such a scenario may be rare, which is why 
we have seen it only once in the UFDs discovered so far. 


3.4. Exotic supernovae 

A possible explanation for high [Fe/H] in a first galaxy 
is a pair-instability supernova (PISN). A single PISN in 
a precursor minihalo could enrich a halo with 10® Mq to 
an average [Fe/H] « —2.5. (Heger et al.|2003l There may 
then be sufficient natural scatter to explain the existence 
of stars at [Fe/H] « —3.5 and —1.5. The suppressed 
[Ba/Fe] and [Sr/Fe] in Segue 1 also supports this ex¬ 
planation. How ever, measurements of a abundances, a s 


well as [Co/Fe] (|Norris et al.||2010al [Frebel et al. 2014|), 
do not match the distinctive signatures expected tor a 
PISN ( Heger fc Woosley]|2002 ). The suggest ion that dif¬ 


ferent elements mix diff erently into the ISM (Sluder et al. 
2015 Ritter et al.|2015 ) may be a partial solution to this. 


however this effect does not explain why all the stars in 
Segue 1 show similar abundance patterns for all a and 
Fe-peak elements. 

Other alternatives to regular core-collapse supernovae 
do not match the abundance patterns of Segue 1. Pul- 
sational pair-instability supernovae eject mostly lighter 
elements, wit h the iron lost as t he star ends its life as 
a black hole (Heger et al.|2003), so are not a solution 
to the observed high [Fe/Hj. Hypernovae can eject more 
iron into the ISM than regular supern ovae, but only by 
a factor of 2 — 3 (Nomoto et al. 20061. They are more 
energetic than regular supernova by an order of mag¬ 
nitude, so will blow out a large proportion of gas and 
metals from a low mass system. Faint supernovae, which 
produce small amo unts of iron compared to carbon have 
been suggested by Cooke & Madau (2014) as an expla¬ 
nation for the formation of CFMP stars m minihalos, as 
their lower energy means they are less efficient at blow¬ 
ing out gas. However, they eject only a small amount of 
iron and would therefore struggle to explain the higher- 
metallicity stars in Segue 1. 

It is therefore unlikely that exotic types of supernovae 
are responsible for the observed abundances of Segue 1. 
However, a pair-instability supernova could produce a 
similar [Fe/H] distribution, so should be considered if 
other systems showing similar MDFs to Segue 1 are dis¬ 
covered. 


4. DISCUSSION 

4.1. The nature and enrichment history of Segue 1 

Segue 1 provides a unique opportunity to investigate 
the early metal build-up in galaxies because it experi¬ 
enced so few chemical enrichment events. In all other 
known galaxies, enrichment from later generations com¬ 
plicates the picture. UFDs were likely some of the early 
building blocks of the Milky Way halo, so by studying 
galaxies like Segue 1, we can help determine the chemi¬ 
cal history of the Milky Way. 

Our results point to a new way of looking at metal 
build-up. When there are few previous enrichment 
events, abundances at a given location can be dominated 
by the effect of the last supernova. [Fe/H] and [C/H] do 
not increase uniformly everywhere in a galaxy, but are 
instead highly inhomogeneous, with some places reaching 
high metallicity before further mixing takes them back 
to close to the average metallicity for the galaxy. This 
“two steps forward, one step back” enrichment means 
that some relatively metal-rich stars form before lower- 
metallicity stars. 

Therefore, while after many supernovae there is a grad¬ 
ual increase in average metallicity over time, this is not 
necessarily the case locally and over shorter periods of 
time. Hence, the oldest stars are not necessarily the 
most metal-poor stars, and more importantly, ancient 
metal-rich stars are expected to exist. None of the sce¬ 
narios that we consider plausible for Segue 1 has the low¬ 
est metallicity stars enriching the gas which forms stars 
with moderate metallicity, which then enriches the gas 
that forms the highest metallicity stars. We find that 
the most likely scenarios have either the lowest metal¬ 
licity [Fe/H] « —3.6 stars enrich the gas that forms 
both the higher metallicity clusters, which share the 
same enrichment history, or that all stars in Segue 1 
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were formed from gas with the same enrichment history. 
Furthermore, given that they require metallicities much 
higher than the likely average of the gas in Segue 1, the 
[Fe/H] « —1.5 stars are more likely to have formed be¬ 
fore the [Fe/H] « —2.4 stars. Gaia will provide accurate 
ages for very metal-poor halo stars and can test for such 
out of order enrichment. If our view of Segue 1 and 
the chemical enrichment of early systems is correct, we 
predict that some of the oldest stars may have higher 
metallicities [Fe/H] > —2.0. 

In systems with few enrichment events, short-term 
variations in environment can result in a large variance 
in [Fe/H]. The first stars in the system can therefore 
be outliers in chemical abundance space. However, if 
star formation is clustered, groups of stars form with 
very similar abundances. The early outliers in abun¬ 
dance space resulting from the varied environment come 
in groups of stars rather than individually. This means 
that rather than each individual system having a small 
number of outlier stars, some systems may have a large 
number (if they contain a cluster or clusters which is 
an outlier in abundance space) and others might have 
none (if the early clusters are not outliers). In this way, 
we explain how 5 of the 7 stars in Segue 1 can have 
A [Fe/H] > 0.8 dex from the mean, whereas most other 
ultra-faint dwarfs (UFDs) have much narrower MDFs. 
The MDFs seen in the majority of UFDs can result either 
from the system having a large number of enrichment 
events, filling in the MDF and hiding distinct clusters, 
or by only having a few clusters with similar [Fe/H]. The 
latter class of systems do not form outlier stars at any 
time. 

For most UFDs, the overall picture is one of early 
environment-dependent enrichment, followed by increas¬ 
ing homogeneity after sufficient enrichment events have 
occurred. In this way, we can account for bot h over¬ 
all homoge neity in the abundances of most stars (Cayrel 


etalj20^^ along with the unusual abundances of others 
([k'rebel (^FNorrisj20I5 Hansen et al.|2015|. The unusual 


stars, which are not necessarily the most metal-poor, are 
those which formed first. 


Our work is in line with that of Ji et al. (2015), who 
suggested that CEMP stars are second-generation stars 
enriched only by massive Pop HI stars. Subsequent gen¬ 
erations quickly wash out this carbon enhancement. This 
fits with Segue 1, where the stars with [Fe/H] « —3.5 are 
CEMP stars, while the higher metallicity stars, which in 
our models are formed from gas enriched by these second- 
generation stars, all have [C/Fe] « 0. Furthermore, it is 
another example of outlier stars being the first to form. 
The yields of Population HI stars appear to be deter¬ 
mined best by studying stars in dwarf galaxies and the 
Milky Way halo which have unusual abundances. 

The discovery of ^ 20 new UFDs in 2015 is likely to 
help answer the question of whether Segue 1 is a rare 
system, or its chemical abundances are common for the 
faintest galaxies. Its MDF is unusual compared to all 
other known UFDs, but nearly all of them are more lu¬ 
minous systems with longer star formation histories. 

4.2. Reticulum 2 

The UFD Reticulum 2 is similar to Segue 1 in its dis¬ 
tance, mean metallicity and luminosity. The system was 
discovered only recently, so at the time of writing only 


iron ( Simon et al. 2015 Walker et al. 2015 Koposov et al. 
2015I)| ~and medium-resolutio n (R ^ 18000) a abun¬ 
dances ( Koposov_et al.|2015b l have been determined. As 
shown in k'igure its MDh' appears similar to the more 
luminous UFDs, contrary to Segue 1. However, this does 
not rule out a similar star formation history to Segue 1. 

In our scenarios above, the [Fe/H] = —2.4 stars in 
Segue 1 formed from gas enriched by the first cluster of 
stars at [Fe/H] = —3.5. This produced a gap in the MDF, 
but equally plausible scenarios may not form such a gap. 
Because of the logarithmic nature of [Fe/H], a first clus¬ 
ter at [Fe/H] = —3.0 would require > 75% of the number 
of supernovae to reach [Fe/H] = —2.5 as a first cluster 
at [Fe/H] = —3.5. It would be much more difficult to 
distinguish between two clusters with means separated 
by 0.5 dex than two clusters separated by 1 dex. 

The scenarios in Section [^ may not always lead to 
the formation of stars withilFe/H] « —1.5 will form. 
Star formation in merging supernova remnants may not 
be common, or the gas whic h returns to the ce nter first 
may not form stars (or, as in Ritter et al. (2015), this gas 
may be carbon-rich and iron-poor). It is therefore possi¬ 
ble that we will hnd many more examples of Ret 2-like 
MDFs than Segue 1, even if both are first galaxies, and 
similar processes occur in both systems. The new galax¬ 
ies discovered by DES provide an excellent opportunity 
to test this. 

There are multiple possible explanations for the origin 
of Ret 2. The first is that it has a similar history to 
Segue 1, with clusters of stars forming at [Fe/H] « —3.2 
and [Fe/H] « —2.5, but no higher metallicity cluster (or 
a small high-metallicity cluster from which we have not 
yet observed stars). Alternatively, the MDF of Ret 2 is 
consistent with a gradual build-up of metallicities as in 
the “Type H” model described in Section 2.1, or at the 
other extreme, a single burst with [Fe/H] « —2.7 and 
cr « 0.3. In Segue 1 we are fortunate that the clusters 
are well-separated, allowing us to distinguish between 
these models. 

Ret 2 is the closest of these new systems, at 30 kpc 
(Koposov et al. 2015a) and is therefore the best candi- 
date for perforrning high-resolution abundance measure¬ 
ments. High-resolution spectroscopy is required to deter¬ 
mine whether Ret 2 passes the tests for a “hrst galaxy” 
as outlined in Frebel & Bromm (2012). A first galaxy 
should not show stars, even at [Fe/HJ > —2, with [a/Fe] 
systematically lower than the halo value of 0.35. Further¬ 
more there should be no signature of s-process or carbon 
enrichment fro m AGB stars. Medium -resolution observa¬ 
tions of Ret 2 (Koposov et al.|2015b) provide a tentative 
suggestion that [a/h'ej does not decline with increased 
[Fe/H], with all 17 observed stars showing [a/Fe] > 0.2, 
close to the mean [a/Fe] = 0.4. High-resolution observa¬ 
tions are required to test this more rigorously, as well as 
to determine the presence or absence of neutron-capture 
elements. 

Segue 1 is highly unusual compared to all other ob¬ 
served galaxies. However, first galaxies may be difficult 
to observe because they are likely to be intrinsically less 
luminous due to their shorter period of star formation, 
so this could be a result of an observational bias. If Ret 2 
proves to be similar to Segue 1, there will be two known 
first galaxies, helping determine which features are typ- 
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ical of first galaxies and which are unique to individual 
systems. 


5. SUMMARY 

We have used the metallicity distribution function of 
Segue 1 to reconstruct its the star formation history and 
to learn about the nature of first and early galaxies. In 
Section]^ we showed that Segue 1 is best explained by 
clusterea star formation. A gradual increase in [Fe/H] 
with time can not reproduce its MDF. There are three 
distinct groups in [Fe/H], with gaps of 1.2 and 0.55 dex. 
All other UFDs observed to date have 90% of their stars 
within 0.5 dex of the mean [Fe/H] for the system and few 
show stars with metallicities as low as [Fe/H] « —3.6 or 
as high as [Fe/H] « —1.5 as in Segue 1. We conclude that 
the unusual MDF is a sign of earliest star formation, in 
line with Segue 1 being a “first galaxy”, meaning that it 
experienced only a single burst of star formation lasting 
< 50 Myr. There are a number of features of Segue 1 that 
are naturally explained by the first galaxy hypothesis: 

1. The gaps in the MDF are caused by the stochas- 
ticity within a single generation of star formation. 
Segue 1 is the only known galaxy which appears not 
to have experienced self-enrichment and therefore 
multiple generations of star formation. Multiple 
generations of star formation are likely to wash out 
stochastic effects, resulting in the MDF being filled 
in near the average metallicity. 

2. The highest metallicity stars remain difficult to ex¬ 
plain, but may be related to the large variation in 
metallicity caused by even a single supernova in a 
low-mass system. First galaxies are more likely to 
have low masses, because they are more suscepti¬ 
ble to effects such as supernova feedback and the 
epoch of reionization, which can quench their star 
formation. 

3. Segue 1 has the lowest intrinsic luminosity of any 
known galaxy. A Kroupa IMF implies that it has 
only formed only « 1500 Mq of stars over its life¬ 
time. This low amount of total star formation is 
likely related to the short duration of star forma¬ 
tion. 

In Section [^ we used hydrodynamical simulations to 
test two poss^le explanations for the MDF of Segue 1. 
One possibility is that the high-metallicity stars formed 
from swept-up gas at the interface between two colliding 


supernova remnants. In this scenario, clusters of stars 
with a difference in [Fe/H] of 2 dex can form within 


10 Myr of each other. Furthermore, as shown by Rit- 


ter et al. (2015) and discussed in Section 2.4, there can 
be large variance in [C/Fe] at early times after a su¬ 
pernova explosion in a low-mass system. In particular, 
the iron-rich ejecta has higher entropy and can overtake 
the carbon-rich ejecta. If this occurs by the time of the 
collision, it provides an explanation for [C/Fe] being sup¬ 
pressed by 1 — 2 dex for the stars at [Fe/H] « —1.5 com¬ 
pared to those at [Fe/H] « —3.5. 

Alternatively, our hydrodynamic simulations show 
that in an Mvir = 10^ Mq dark matter halo with 10® Mq 
of gas pre-enriched to [Fe/H] = —2.5, the first gas that re¬ 
turns to the center of the system after a supernova explo¬ 
sion will be enriched to [Fe/H] « —1.5 and be sufficiently 
dense that star formation is plausible. We also briefly 
explored the possibility of more exotic types of super¬ 
novae such as PISNs to reproduce the MDF of Segue 1, 
but these could not explain the abundances of other ele¬ 
ments. 

At the time of writing ^ 17 new UFD candida t es have 
been detected in DES (Koposov et al.||2015a| |Drlica-j 


Wagner et al. 2015). Oiie of these. Reticulum 2, which 
has been confirmed as a galaxy, is at a similar distance 
as Segue 1 and has a similar mean metallicity. Medium- 
resolution observations of [a/Fe] suggest it as a tentative 
candidate for a first galaxy. High-resolution abundance 
measurements are required to determine whether this is 
indeed the case. Segue 1 is the only known first galaxy 
and it is therefore currently impossible to say whether 
its features are typical or unusual. The newly-discovered 
systems provide an exciting opportunity to test this. In 
the long-term, LSST will be capable of detecting much 
fainter objects, while G-CLEF on the GMT will enable 
high-resolution abundance measurements of fainter stars. 
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